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A BAYESIAN DISCRIMINATOR FOR RAPIDLY DETECTING ARRHYTBQMIAS 

1. FIELD OF THE INVENTION 

The invention is directed to the generation and analysis of data with a multiple-index 
Bayesian discriminator. More spedfically, the invention is directed to methods, systems, and 
5 devices for drtecting and treating arrhythmias and heart diseases. 

2. BACKGROUND OF THE INVENTION 
2.1 Arrhythmias 

Arrhythmias are caused by a disruption of the normal electrical conduction system of 
the heart, causing abnormal heart rhythms. Normally, the four chambers of the heart (2 atria 

.0 and 2 ventricles) contract in a very specific, coordinated manner. The signal to contract is an 
electrical impulse that begins in the "sinoatrial node" (the SA node), which is the body's 
natural pacemaker. The signal then travels through the two atria and stimulates them to 
contract. The signal passes through the "atrioventricular note" node (the AY node), and 
finally travels through the ventricles and stimiilates them to contract. Problems can occur 

15 anywhere along the electrical conduction system, causing various arrhythmias. There can be 
a problem in the heart muscle itself, causing it to respond differently to the signal, or causing 
the ventricles to contract out of step with the normal conduction system. Other causes of 
ardiythmias include abnormal rhythmicity of the body's natural pacemaker, a shift of the 
pacemak^ fi-om SA node to other parts, blocks at different transmission points, abnormal 

20 pathways of impulse conduction, and spontaneous gmeral of abnormal impulses due to 
ischemia (low flow to coronary arteries), hypoxia (low oxygen), ANS imbalance, lactic 
acidosis, electrolyte abnormality, drug toxicity, and hemodynamic abnormalities. 

Atrial fibrillation (AF) is the most common form of supraventricular arrhythmia and is 
associated with a considerable risk of morbidity and mortality. (Benjamin EJ, et al., 1998 

25 C/roito/ow 98:946-952; Ryder KM, etaL, 1999 Am J Cardiol 84:1311R-138R; ChughSS, et 
al. 2001 J Am Coll Cardiol 37:371-377). As many as 2 million Americans are Uving with 
atrial fibrillation according to the American Heart Association Theoretical analyses and 
high-density mapping studies have suggested that the most common mechanism of AF is the 
presence of multiple wave fronts or "wavelets" circulating irregularly throughout the atrial 
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tissue. (Moe GK, et al., 1964 Am Heart I 67:2961-2967; Allessie MA, et al., "Experimental 
Evaluation of Moe s Multiple Wavelet Hypothesis of Atrial Fibrillation" in Zipes EP, Jalife J, 
eds. Cardiac Electrophysiology ondArrhyhtntias, Orlando, Fla: Grune & Stratton, Inc., 1985; 
pp 265-275; Konings KTS, et al., 1994 Circulation 89:1665-1680). 

5 Various studies have employed lime domain, frequency domain, or time-frequency 

analysis to differentiate fibrillatory from non-fibrillatory rhythms. However, most of these, 
essentially single-index methods, techniques suffer from limitations such as rather long 
process time, lack of robustness to noise or far field events, poor performance in 
discriminating atrial flutter (AFL) from sums rhythm (SR), and relatively low sensitivity and 

10 specificity. 

These also limit improvements in pacemakers and other devices. For instance, in 
dual chamber pacemakers, accurate AF detection is critically important to avoid rapid 
ventricular pacing by activating automatic mode switching. In an implantable cardioverter 
defibrillator, accurate recognition of AF can avoid false discharges. Furthermore, the recent 

15 development of automatic implantable atrial defibrillators has created a critical need for 
speedy and reliable disOTmination of AF from other types of intra-atrial electrograms. (Lau 
CP, et al., 1997 Pacing Clin Electrophysiol 20:220-5; Wellens HI, et al., 1998 Circulation 
98:165 1-1656; FriedmanPA, et al., 2001 Circulation. 104:1023-1028). 

Proposed techniques for detecting AF can be conveniently divided into about four 

20 categories such as (1) methods based on time domain features (See Botteron GW, et al., 1996 
Circulation 93:513-518; Botteron GW, et aL, 1995 IEEE Trans. BME 42:579-586; Tse HF, et 
al., 1999 Circulation 99:1446-1451; Sih HJ, et dL, 1999 IEEE Trans. BME 46:440-450; 
Swerdlow CD, et aL, 2000 Circulation 101:878-885; Thakor NV, et oL, 1990 IEEE Trans. 
BME 37:837-843; Chen SW, et aL, 1995 J Electrocardiol 828:162; Chen SW, et al., 1996 

25 TREK Trans. BME 43:1120-1125); (2) methods employing frequency domain properties (See 
RopeUa KM, et aL, 1989 Circulation 80:112-119; Bollmann A, et aL, 1998 Am J Cardiol 
81:1439-1445; Chen SW, et al., 1996 ICASSP P(J 3:1775-1778; Chen SW, 2000 IEEE Trans. 
BME 47:13 17-1327); (3) techniques making use of tune and frequency analysis (See Slocum J, 
et al Computer discrimination of atrial fibrillation and regular atrial rhythms from intra- 

30 atrial electrograms. Pacing Clin Electrophysiol 1988;11:610-621; Lovett EG, et al, 1997 
Ann BME 25:975-984); and (4) miscellaneous (See Zhang XS, et al, 1999 IEEE Trans. BME 
46:548-555). 
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Botteron and Smith developed an algorithm based on the crosscorrelation of two pre- 
processed bipolar intra-atrial signals of which an active space constant was extracted (1996 
Circulation 93:513-518; 1995 IEEE Trans. BME 42:579-586). Tse et aL depicted a two- 
phase AF detection method that directly processed the time domain signals (1999 Circulation 

5 99:1446-1451). More recently, Sih et aL proposed an approach employing the mean square 
error in the linear prediction between two unipolar epicardial electrograms (1999 IEEE Trans. 
BME 46:440- 450). Swerdlow et oL used a technique that combmed the median cycle length 
and an atrial tachyarrhythmias evidence counter that used the number of sensed atrial 
electrograms in consecutive RR intervals (2000 Circulation 101:878-885). Chen et aL 

10 proposed a modified sequential algorithm based technique (1995 J ElectrocardioL 828:162; 
1996 IEEE Trans, BME 43:1120-1125). Instead of measuring the rate, they employed 
blanking variability to measure the temporal irregularity with improved detection accuracy. 

In addition to the time domain measures mentioned above, there are methods rooted in 
spectral analysis including coherence spectrum method and frequency analysis using the 

15 surface electrocardiogram. (See Ropella KM, etal, 1989 Circulation 80:1 12-1 19; BoUmann 
A, et al,, 1998 Am J Cardiol 81:1439-1445). Chen et aL disclosed a two-stage arrhythmia 
discrimination method using a damped e^onential modeling algorithm which gives higher 
frequency resolution than simple Fast Fourier transform methods (1996 ICASSP 96 3:1775- 
1778; 2000 TKFK Trans. BME 47:1317-1327). Similarly, Slocum et al designed an 

20 algorithm that took into account both the morphological information (atrial rate and amplitude 
probability function) and frequencgr domain features (power spectrum analysis) (1988 Pacing 
Clin ElectrophysioL 1 1 :6 10-621). In addition, Lovett and Ropella disclosed analysis of atrial 
riiythms via a time-frequency distribution of coherence (1997 Ann BME 25:975-984). From 
the viewpoint of dynamical systems, Zhang et aL proposed a complexity-based approach for 

25 discrimination of ventricular tachycardias and fibrillation (1999 IEEE Trans. BME 46:548- 
555), a method having a few advantages over the conventional detection techniques (Chen SW, 
2000 IEEE Trans. BME 47:1317-1327). However, these methods need rather long episode 
(>5s) to get satisfactory performance. 

While techniques using single-index calculation are useful in the detection of 

30 arrhythmias, there is a continued need to find more accurate and rapid detection modalities 
and approaches to diagnose arrhythmias. 



2.2 Statistical Analysis 
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2.2.1 Bayes Decision Rule 

A Bayesian theorem describes the relationship that exists between simple and 
conditional probabilities. The Bayes decision theory assumes that the decision problem 
(whether an observed episode belongs to one class or another) is posed in probabilistic terms, 
S and that all of the relevant probabilities are known. For Instance, P(Wj) is denoted to be the 
prior probability that a certain episode should belong to W|, i.e., P(w|) is the probability that 
an episode is of class i even before it is observed. The symbol p(v | denotes the class 
conditional probability of observing feature vector y given the fact v is of class Wi is known. 
In other words, p(v 1 wj is a probability density function of non-negative value and can be 
10 estimated by the trainmg data set. P(Wi | v ) is called the posterior probability which can be 
calculated by p(v | Wi) and P(Wi) according to the Bayes' rule. P(Wi | v) is the probability 
(between 0 and 1) that an object is of class Wj given it is observed as v .If the cost of a correct 
decision is 0, and the cost of a wrong decision is 1, then, the Bayes Decision Rule can be 
applied as: Decide Wi P( w^ | v ) > P( Wj| v ) for all 

IS 2.2.2 Sensitivity, Specificity, and Accuracy 

Sensitivity and specificity together describe the accuracy of a test. When a large 
number of positive and negative samples are tested, sensitivity determines the percentage of 
false-negative results, and specificity determines the percentage of false-positive results. For 
example, a specificity of 99% means that 1% of those without AF will test felse-positive for 
20 ejdiibiting AF. A sensitivity of 99%, on the other hand, means that 1% of those with AF will 
test false-negative, /.e., as not cTdiibiting AF. 

3. SUMMARY OF THE INVENTION 

Atrial fibrillation (AF) is the most common arrhythmia (abnormal heart beat) with a 
considerable risk of stroke and mortality. Atrial flutter (AFL) is another type of abnormal 
25 heart beat that also occur frequentiy in those patients with AF. Accurate and rapid detection 
of these rhythms is critically important to avoid rapid ventricular pacing by activating 
automatic mode switching and false shock discharges from implantable device (pacemaker 
and defibrillator). The detection of these abnormal rhythms by implantable devices require 
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the use of intra-a^al electrograms recx)rded from the atria. Since the treatments of AF & AFL 
are clinically completely different, it is of rather urgent need that an algorithm is written to 
distmguish these three types of heart signals by the device. To train up our system of 
detection, several hundred episodes of intra-cardiac signals (called the closed data set CDS) 

5 were recorded by a computer. Five feature parameters were evaluated for each episode or 
window, and a discriminator is obtained to decide which class of signals (AF, AFL or smus 
rhythm (SR.), normal heart beat) does this episode belong to according to the mathematical 
method specified below. Experienced physicians make also a dedsion for each episode 
independently. The two results are then compared. The performance of this algorithm as 

10 specified by the specificity, accuracy and sensitivity. After checking these three to be 
satisfactory (>97%), the statistical averages of the five feature parameters are calculated and 
the system is ready to use. 

Our new algorithm allows decision to be made based on a window of intracardiac 
signal of interval 1-2 seconds only and the computer calculation time is shorter than 0.25s. 

15 Any new set of episode is added to the CDS. Since the performance of tiie algorithm 
improves as the size of CDS is increased, this algorithm gets "smarter" as more cases are 
tested. 

To check the performance, we have used several hundred episodes of rhythms called 
open data set (ODS, different from CDS) and have found that our methodology works well. 

20 We impose noise to the ODS and found that the algorithm has very good anti-noise property. 

Note that we have used five feature parameters based on the physical interpretation 
specified later. This number five can be extended to higher number for better result, when 
we find other interpretations or when we treat signals other than AF, AFL & SR. Moreover, 
based on information from one or few windows ('^ 1 second) of signals, the calculation time 

25 has to be very short (preferably < Is) so that the implantable device can use the information 
and make decision on the type of treatment on line. Our invention marks the basis of 
produdng software to be attached to machmes associated with intracardiac signal detection. 

4. BRIEF DESCRIPTION OF FIGURES 

FIGS. lA-C illustrate the various feature extraction for episodes of SR (FIG. lA), AF 
30 (FIG. IB), and AFL (FIG. IC) including (a) raw episode; (b) output after manipulations 1 to 
3; (c) auto-corrdation coeflScients; and (d) rectified version read for feature extraction. 
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FIG. 2 shows a flowchart of the steps mvolved m the traming and discrimination 
procedure. Block arrows indicate the training process and solid arrows indicate the detection 
procedure. 

FIG. 3 shows the comparison of vahies of features for open and close data sets. 
5 White, black, and shadowed bars represent SR, AFL, and AF, respectively. Each of the five 
features are significantly different between AF, AFL, and SR for both open and close data sets. 
There are no significant differences in the values of each of the five features for AF, AFL, and 
SR between close and open data set 

FIG. 4 shows the performance (e.g., sensitivity, specificiQr, and accuracy) achieved 
1 0 according to the number of features used. 

BIG. 5 shows the relationship between the performance (e.g., sensitivity, specificity, 
and accuracy) of the disclosed dismminator and the signal-to-noise ratio (SNR). 

5. DETAILED DESCRIPTION OF THE INVENTION 

The present invention generally relates to methods, systems, and devices for 

15 detecting and treating arrhythmias and heart diseases. Atrial tachyarrhythmias are detected 
in a subject using a multiple-index Bayesian discriminator. The method for detection 
comprises the steps of obtaining an open-test data set of bipolar intra-atrial signals from the 
subject of interest and using a computer or computers to analyze the open-test data set. 
Furthermore, the method for detection generates a result in accordance with a set of estimated 

20 conditional probabilities from a trmiing data set based on the multiple-index Bayesian 
discriminator. The use of a computer, or a computing device system in practicing the method 
is illustrative and includes any computer executable processing device. Similarly, the 
method is suitable for detecting various conditions such as sinus rhythm, atrial flutter, atrial 
fibrillation, or any type of arrhythmias, heart diseases, or physiological conditions. In 

25 general, the open-test data set may comprise any type of electophysiological information (e.g,, 
ECG, EEG, and EKG) obtained from the subject of interest althougfi ECG data is employed in 
the preferred embodiment. 

In another embodiment, the method for detection further comprises the steps of 
selecting a phirality of features of intra-atrial dectrograms and a ^e of output, inputting a 

30 close-test data set of bipolar intra-atrial signals for training, and estimating tiie set of 
conditional probabilities for the plurality of features and the type of output m accordance with 
a multiple-index Bayesian discriminator from the close-test data set. Of course, the method 
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described herewith is applicable to any type of electophysiological information (e.g., ECG, 
EEG, and EKG) obtained from the subject of interest. 

In another embodiment, the method for detection fiirther comprises the step of 
selecting additional features for estimating conditional probabilities. The phirality of features 
5 of intra-atrial electrograms may be selected from the non-exhaustive illustrative list 
comprising regularity, rate, energy distribution, percent tune of quiet interval, and number of 
baseline reaching. For instance, the phiralily of features may also be selected from those 
parameters disclosed in previous studies such as cross-correlation of two pre-processed 
biopolar intra-atrial signals (Botteron GW and Smith JM, 1995 IEEE Trans. EME 42:579-586; 

10 Botteron GW and Smith JM, 1996 Circulation 93:513-518), time (Tse HF, et al., 1999 
Circulation 99:1446-1451; Thakor NV, et al., 1990 IEEE Trans, BME 37:837-843), mean 
square error in the linear prediction between two unipolar epicardial electrograms (Sih HJ, et 
al., 1999 IEEE Trans, BME 46:440- 450), median cycle length in conjunction with the number 
of sensed atrial electrograms in consecutive RR mtervals (Swerdlow CD, et aL, 2000 

15 Circulation 101:878-885), temporal irregularity (Chen SW, et al, 1995 J Electrocardiol. 
828:162; Chen SW, et al., 1996 IEEE Trans. BME 43:1120-1125), and frequency (Ropella 
KM, et al. 1989 Circulation 80:112-119; BoUmann A, et al. 1998 Am J Cardiol 81:1439-1445; 
Chen SW, et al., 1996 ICASSP 96 3:1775-1778; Chen SW, 2000 IEEE Trans. BME 47:1317- 
1327). 

20 In another embodiment, the method for detection iiirther conq)rises the step of 

modifying at least one estimated conditional probabilities from the set of estimated conditional 
probabilities. Preferably, the open-test data set and the results obtained from analysis of the 
open-test data set are incorporated into to the closed-test data set in an iterative maimer. The 
set of estimated conditional probabilities is continuously modified as more data set is inputted. 

25 Tfaus^ performance of the method can be contiimously modified or improved, increasing 
the specificity, sensitivity, and accuracy of the result. 

In another embodiment, the method for detection fiulher comprises the step of 
differentiating between the types of aithythmias or heart diseases in the subject of interest 
To this end, a sufficient nmnber of features of intra-atrial electrograms are used so the method 

30 for detection displays an overall sensitivity of at least 90%, preferably 95%, more preferably 
98%, and most preferably 99%, an overall specificity of at least 90%, preferably 95%, more 
preferably 98%, and most preferably 99%, and an ovearall accuracy of at least 90%, preferably 
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95%, more preferably 98%, and most preferably 99%. An illustrative non-exhaustive Ust of 
arrhythmias detected by the disclosed method includes sinus rhythm, atrial flutter, atrial 
fibrillation, atrial tachyarrhythmias, tachycardia, bradycardia, supraventricular airhythmias, 
premature atrial contractions (PACs), paroxysmal supraventricular tachycardia (PSVT), 
5 accessory pathway mediated tachycardias, atrial tachycardia, ventricular arrhythmias, 
premature ventricular contractions (PVCs), ventricular tachycardia, ventricular fibriUation, 
bradyanhythmias, sinus node dysfunction, and heart block. 

In another embodiment, the method for detection shows robust anti-noise performance 
in differentiating between atrial fibrillation (AF), atrial flutter (AFL), and sinus rhythm (SR). 

10 The overall sensitivity, specificity, and accuracy of a method for detection is similar at 
diflferent signal-to-noise ratio (SNR) above 10 dB. The overall sensitivity of the method for 
detection is at least 90%, preferably 95%, more preferably 98%, and most preferably 99% 
when the SNR is greater than 10 dB. Similarly, the overall specificity of the metiiod for 
detection is at least 90%, preferably 95%, more preferably 98%, and most preferably 99% and 

15 the overall accuracy of the method for detection is at least 60%, 65%, 70%, 75%, 80%, 85%, 
90%, or 95% when the SNR is greater than 10 dB. 

In another embodiment, the method fiirther comprises the step of providing a treatment 
in response to detecting a particular conditioa Such treatment options include, but are not 
limited to, medications, cardioversion, pacemakers, implantable cardioverter-defibrillators, 

20 surgery, or radiofirecpiency catheter ablation of the anhythmia focus. In particular, an 
implanted device that can adjust its stimulation in response to rapidly detecting a particular 
arrytiimia. Such rapid detection is cabled in less than five seconds, more preferably in less 
than 4 seconds, even more prrferably less than 3 seconds and most preferably less than 2 
seconds mcluding at least one of 1.9 sees., 1.8 sees,, 1.7 sees., 1.6 sees., 1.5 sees., 1.4 sees., 

25 1.3 sees., 1.2 sees, 1.1 sees., 1.0 sees., 0.9 sees., 0.8 sees., 0.7 sees., 0.6 sees., 0.5 sees., 0.4 
sees., 0.3 sees., 0.2 sees, and 0.1 sees. 

In another embodiment, a device detects arrhythmias in a subject of interest. The 
device for detection comprises a module for collecting an open-test data set of bipolar intm- 
atrial signals fi-om the subject of interest and a computer or a system of computer devices for 

30 analyzmg the open-test data set Furthermore, the device for detection comprises a screen or 
similar device that can display the results ia accordance with a set of estimated conditional 
probabilities. The open-test data set can be collected in any tangible or intangible database or 
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Storage means. The module need not be a separate or discrete unit; it can be a program, a 
processor, a sub-component, etc. Further, the analysis could be carried out by any computer 
executable processing device and not just a computer. Similarly, the device could be used to 
detect sinus rhythm, atrial flutter, atrial fibrillation, or any type of arrhythmias, heart diseases, 
5 or physiological conditions. 

In another embodiment, the device for detection further comprises a module, wherein 
the module selects a plurality of features of intra-atrial electrograms and a type of output, 
inputs a dose-test data set of bipolar intra-atrial signals for training, and estimates the set of 
conditional probabilities for tiie phirality of features and the type of output in accordance with 

10 a multiple-index Bayesian discriminator from the dose-test data set. The device for 
detection furth^ comprises a tUrd module, whearein the module selects additional features for 
estimating conditional probabilities. Possible features of intra-atrial electrograms for 
analysis include the features in tiie group consisting of regularity, rate, energy distribution, 
percent time of quiet interval, and number of baseline reacMng, cross-correlation of two pre- 

15 processed biopolar intra-atrial signals, time, mean square error in the linear prediction between 
two unipolar epicardial electrograms, median cycle length in conjunction with the number of 
sensed atrial electrograms in consecutive RR intervals, temporal irregularity, and frequency. 

A module may perform all or a sub-combination of steps, collecting data set, 
analyzing data set, providing an analysis, selecting a plurality of features, selecting a type of 

20 output, estimatmg a set of conditional probabilities, and displaying the intermittent and/or final 
results. Further, the analysis could be carried out by any compute executable processing 
device or devices. Furthermore, the module may inchide facility for modification of an 
estimated conditional probabilities from the set of estimated conditional probabilities. In 
order to so modify any conditional probability, preferably, the open-test data set and the 

25 results obtained from analysis of the open-test data set are added to the closed-test data set m 
an iterative manner. The set of estimated conditional probabilities is continuously updated as 
more data set is inputted. Thus, the performance of the method is continuously modified or 
improved, /.e,, increasing the specificity, sensitivity, and accuracy of the result Of course, 
more than one estimated conditional probabilities may be improved upon in like manner. 

30 In another embodiment, a device for detection further comprises a module, wherein the 

module differentiates between the types of arrhythmias or heart diseases in the subject of 
interest. In a preferred embodiment, the module uses a sufficient number of features of intra- 
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atrial electrograms so the device for detection displays an overall sensitivity of at least 90%, 
preferably 95%, more preferably 98%, and most preferably 99%, an overall specificity of at 
least 90%, preferably 95%, more preferably 98%, and most preferably 99%, and an overall 
accuracy of at least 90%, preferably 95%, more preferably 98%, and most preferably 99%. 

5 The different types of arrhythmias include, without limitation, sinus rhythm, atrial flutter, 
atrial fibrillation, atrial tachyarrhythmias, tachycardia, bradycardia, supraventricular 
anhythmias, premature atrial contractions (PACs), paroxysmal supraventricular tachycardia 
(PSVT), accessory pathway mediated tachycardias, atrial tachycardia, ventricular arrhythmias, 
premature ventricular contractions (PVCs), ventricular tachycardia, ventricular fibrillation, 

10 bradyarrhythmias^ sinus node dysfunction, and heart block. 

In yet another embodiment, a device for detection further comprises a member that 
provides a modulating effect on heartbeats corresponding to the result For instance, the 
member can deliver an electrical signal or input to the chest wall that synchronizes the heart 
and allows the normal rhythm to restart (as in a electrical cardioversion). Or, the member 

15 can send small electrical impulses to the heart muscle to maintain a suitable heart rate (like a 
pacemaker), deliver energy to the heart muscle to cause the heart to beat in a normal rhythm 
(like an implantable cardioverter-defibrillator), and even direct applying or delivering of high 
radio-firequency energy through a special catheter to small areas of tissues that cause abnormal 
heart ihytiims (as in radiofirequenc^ catheter ablation). Moreover, this description of the 

20 memb^ is illustrative rather than limiting. For instance, different types and combinations of 
pacemakers and implantable cardioverter-defibrillators can be directly incorporated into the 
device. Additional technology for modulating (i.e., increases, decreases, stabilizes) heart 
rhythms can be incorporated into the device without limitation to respond to the detection of a 
particular anhythmia. Such technology can include pharmaceutical, biological, chemical, 

25 physiological, electrical, anatomical, and molecular (Le., antibodies, anti-antibodies, fusion 
protdns, polypeptides, firagments, homologues, derivatives, and analogues thereoQ 
possibilities. 

The subjects to which the methods, systems, and devices for detection and treatment of 
the present invention are applicable may be to any mammalian or vertebrate species, which 
30 include, but are not limited to, cows, horses, sheep, pigs, fowl (e.g., chickens), goats, cats, 
dogs, hamsters, mice, rats, monk^s, rabbits, chimpanzees, and humans. In a preferred 
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embodiment, the subject is a human. Additional teachings are clarified with the aid of details 
in an example study below. 

5.1 EXAMPLES 

5.1.1 Data Acquisition 

5 Bipolar mtra-atiial electrograms at high anterolateral right atrium (with a 1 cm inter- 

electrode distance) from 20 patients in AF, AFL and SR were amplified and recorded 
(CardioLab 4.11, Pruka Engineering, Inc.) during electiophysiologicai procedures. The 
patients were presented to the electrophysiology laboratory for internal cardioversion of AF, 
electrophysiology study and/or radiofrequency ablation procedure for their underlying 

10 arrhythmias. Up to 220 seconds (mean: 190 ± 20 seconds; range: 180 to 220 seconds) of 
simultaneous imfiltered (band pass 0.04-SOOO hertz) recording from each patient were 
digitized at 1000 hertz. The data was then split into 1 (AF & AFL) or 2 seconds (SR) 
segments for analysis so that at least two atrial events were recorded during SR. In order to 
geaerate an unbiased data set, nearly the same numbers of episodes were randomly collected 

15 &Qm each patient. Computer processing was performed using a Matlab 5.3 computer 
program (The Mathwork, Inc.). 

The example study consisted of 20 patients (17 men and 3 women, mean age 55±16 
years, ± SD). Their mean left ventricular ejection fraction was 56±10%, and their mean left 
atrial diameter was 4.6±1.7cm as measured by echocardiography. Their clinical 

20 characteristics are summarized in TABLE 1. 



TABLE 1, Patients Characteristics 



Patient 


Age 


Sex 


Diagnosis 


Medications 


Rhythm 

recorded 


Procedure 


1 


50 


M 


HT.AF 


CCB, Amiodarone 


SR, ST, AF 


InteraalCV 


2 


55 


M 


Lone AF 


Amiodarone 


SR,AF 


Internal CV 


3 


68 


M 


HT.AF 


BB, Amiodarone 


SR,AF 


IntemaiCV 


4 


55 


M 


HT.AF 


ACEI^ Amiodarone 


SR,AF 


IntonalCV 


5 


60 


M 


Lone AF 


CCB 


SR,ST,AF 


IntemaiCV 


6 


53 


M 


HT.AF 


ACEI, Amiodarone 


SR,AF 


IntemaiCV 
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7 


72 


F 


Congestive heart 
failure, AF 


ACEl, Digojdn, 
Amiodarone 


SR,AF 


Internal CV 


8 


48 


M 


Pericarditis, AFL 


Sotalol 


SR, Typical 
AFL 


EP/RF 


9 


53 


M 


Coronary artery 
disease, AFL 


BB,CCB 


SR,ST, 
Typical AFL 


EP/RF 


10 


70 


M 


Lone AFL 


None 


SR, 
Atypical AFL 


EP/RF 


11 


66 


M 


HT.AF 


CCB, Amiodarone 


SR,AF 


Internal CV 


12 


67 


M 


HT.AFL 


CCB 


SR, Typical 
AFL 


EP/RF 


13 


64 


M 


Lone AFL 


CCB 


SR, Atypical 
AFL 


EP/RF 


14 


40 


M 


Lone AF 


Amiodarone 


AR,AF 


Internal CV 


15 


66 


M 


HT.AF 


CCB, Amiodarone 


SR,ST 


Internal CV 


16 


50 


M 


AVNRT 


None 


SR,ST 


EP/RF 


17 


56 


F 


WPW 


None 


SR,ST 


EP/RF 


18 


14 


F 


WPW 


None 


SR,ST 


EP/RF 


19 


21 


M 


WPW 


None 


SR,ST 


EP/RF 


20 


70 


M 


AVNRT 


None 


SR,ST 


EP/RF 



Abbreviations: AF, atrial fibrillation; AFL, atrial flutter, AVNRT, atrioventricular nodal 
reentrjr, BB, beta-blocker; CAD, coronary artery disease; CCB, calcium channel blocker, EP, 
electrophysiology study; HT, hypertension; RF, radiofrequency ablation; SR, sinus rhythm; 
5 ST, sinus tachycardia; WPW WolfF-Parkinson-White syndrome. 

A total of 364 bipolar recording were collected from these patients. All rhythm 
episodes have been assessed blindly and classified mto AF, AFL or SR by 2 experienced 
electrophysiologists. Of these recording, 156 episodes were AF, 88 episodes were AFL 
(mean atrial cycle length 320±40 ms, range 290-345 ms), and 120 episodes were SR, 
10 mcluding 50 episodes of sinus tachycardia during isoprenaline infiision (mean sinus cycle 
length 535i30 ms, range 505-570 ms). Eadi patient contributed nearly the same number of 
qjisodes to the data set (18-22 episodes per patient). We randomly selected 219 (60%) and 
145 (40%) riiythms as close-test data set and openrtest data set, respectively. 
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5.1.2 Signal Manipulation 

Before extracting the features of the signal, each rhythm episode was processed wth 
the foUowing manipulations: (1) third-order Butterworth bandpass filtering (40-250 Hz), (2) 
absohite vahiing, (3) low pass filtering (0-20 Hz), (4) autocorrelation, and (5) rectification 
5 (FIGURE 1). Steps I to 3 ou^ut a flattened signal proportional to the high frequency 
energy contamed in the mput episode. (Botteron GW and Smith JM, 1996 Circulation 
93:513-518; Botteron GW and Smith TM, 1995 IEEE Trans. BME 42:579-586). The 
autocorrelation process avoids drastic fluctuation of the ampUtude of atrial electrograms with 
time. (Oppendheim AV, Schafer RW. In: Discrete-Time Signal Processing, Chapter 11, 
10 Prentice-Halllntemational, Inc., 1989:742-756. Krauss TP, Shore L, Little JN. In: Signal 
Processing Toolbox User's Guide, Chapter 1, The Matii Works Inc., 1994:61-63). Finally^ 
tiie rectification process removes all tiie negative parts of tiie processed signal to fecilitate the 
mathematical treatment during feature extraction. 

5.1.3 Feature Extraction Procedure 

15 Five relevant feature parameters were extracted fi-om the final processed signal by a 

feature extraction procedure (FIGURE 1). The first feature (f^ is defined as the first peak, 
occurrmg at time (t), which is positively related to tiie regularity of tiie mput The second 
feature (Q is defined as f2= t/1000, and is proportional to tiie input's atrial rate. The tiiird 
feature (fj) is defined as tiie percentage of energy contamed m flie two time bands (Ei+Ej^), 

20 where El, is the energy witiun 0 to 100 ms, E, is tiie energy witiiin 500 ms to 1000ms, and E is 
tiie total energy witiun 0 to 1000 ms. The typical suras rate is measured at 60-120 beats per 
mmutes, flie corresponding peak to peak interval is 500-1000 ms. In SR, tiie energy is 
mainly distributed m the aforementioned two time bands. Therefore, feature 4 is helpfiil to 
distinguish SR signals from tiie otiier two classes of rhytiua (AF or AEL) since the vahie of fj 

25 is very close to one for SR and smaller for AF or AFL. The fourth feature (Q measures tiie 
percent time interval corresponding to zero amplitude signal (percent quiet mterval) and is 
calculated by tiie sum of time mtervals witii zero vahie over tiie total duration of rectified auto- 
conetetion fimctioa The fifth feature (fs) measures tiie number of components tiiat reachmg 
tiie baselme in 1 second (baseline reaching). Botii features f* and ^ reflect tiie chaotic extent 

30 or randomness of tiie input signals and flierefor, are supposed to be sensitive to fibrillatory 
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rhythm (AF). The entire group of parameters ^, fg, and ^ form a vector in five 
dimensions, which can only be determined if all the values of these 5 variables are knowa 

FIGURE 4 shows respectively the sensitivity, specificity and accuracy of rhythm 
detection versus the increase of features. With the number of feature(s) used increase fi"om 1 to 
5 5, the performance increases significantly Q)<0.01) firom aroxmd 80% to above 95. This result 
also indicates the advantage of multi-feature detection over single-feature detectioa 

The results of S extracted features for the close and open data set are presented in 
FIGURE 3. The values of each of J features were significantly different between AF, AFL 
and SR for both close and open data set. However, there are also significant overlaps 
10 between the values among the three types of rhythm for each feature. There were no 
significant difBa-ences m the values of each of 5 features for AF, AFL and SR between close 
and open data set, suggesting the two data sets were very similar. 

5.1,4 Training Process 

Sixty percent of the collected rhythm episodes were randomly selected as the closed- 
15 test training data set of the new discriminator. The values of f^, fz, fa, and f^, and the 
corresponding feature vector for the three classes of rhythm signals (SF, AF, and AFL) were 
obtained. The distribution of each of the five features has been found to follow 
approximately the normal distribution, therefore, the corresponding feature vectors of each 
class of rhythm also satisfy approximately a 5-dimensional normal distribution. Similar to 
20 the one-variate normal distribution, the multi-vanate normal distribution is also determined by 
two parameters - mean and the so-called Covariant Matrix, both of which could be estimated 
via the feature vectors of the training data set (close-test data set). The mean and the 
Covariant Matrix are both necessary for the discrimination procedure depicted in the following 
section 6.1,5. 

25 The objective of training process is to estimate the prior probabilily P(Wj) and the class 

distribution p(v | Wj) . These two items are necessary for calculating the posterior probability 
p(Wil v) which is critical for the discrimination procedure. In practice, P(Wi) can be 
approximated by /2,/Z/=i W/, where is tiie total episode number of the i^ class. P(Wi | v)can 
be calculated by p( v | Wj) P(Wi) according to Bayes' rule. Assume that p(v | wj is normal, 

30 the foUov^g equation (1) is obtained: 
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'^^^'^^^^^(^'^^ 

where /i= E [v ] is the mean of v, and 
5 E=E[(v- /i)(v - y5)*] is Ihecovariant matrix generated by Ihe vector (v /i);t denotes 
transpose and -1 denotes inverse of a matrix. 

5.1.5 Discrimination Procedure 

In order to optimize detection performance, a multi-variate Bayes decision theory is 
used. (See section 2.2.1). Using the Bayes Theorem, the posterior probability, which is the 
10 chance that a feature vector of any episode should belong to any of the three classes of rhythm, 

is calculated. Then, a so-called "discrimination function, g (v)" or a class of rhythm in 
general based on Bayes decision theory, is generated. For each rhythm episode, the values of 
the three discaimination functions gg^ (v), g^p 00, gAPLCO* which correspond to the probabilities 
of the episode belonging to SR, AF, and AFI^ respectively, are evaluated. The final decision 
IS for each rhythm episode is simply determined by which of absolute value of the above three is 
the largest (FIGURE 2). The detailed mathematical treatment leading to the representation 
of the discrimination function is discussed below. 

Theoretically, the detection process is to calculate the posterior probabilities 
P(Wi| v)= p(v|Wi)P(Wi) of all 3 classes given one unknown episode. However, because 
20 normal distribution has e?q)onential terms, which is time-consuming to calculate, for 
computation e£&ciency, the logarithm on both side of the above equation is taken: 

gi(v = logP V I Wi) + logP(Wi) (2) 
Then, equation (1) is substituted into equation (2), obtaining a convoiient form for the 
"discrimination function" &( v ): 
25 gi(v)= v^Wi v + w> i7+Wi„ (3) 

whereWi = -V42:;^ (4) 
w = V A (5) 
Wio = -i^2A*2:i-^A-VilogeISi| + lo&P(w^ (6) 
After calculating the three values of &v (i = 1,2,3), the i value corresponding to the ma x im um 
30 g| is chosen according to the Bayes decision rule. 
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5.1.6 Anti-noise Performance 

Sometimes the intracardiac signals may be comipted by noises introduced by external 
electromagnetic interference and myopotential sensing. It is important for the method to be 
robust \^en processing noisy episodes. As shoTvn in this study, the SNR. has significant 

5 effect on the performance of the disclosed discriminator. A deorease in SNR reduces the 
sensitivity for detection of regular rhythms, such as SR and AFL. This phenomenon is due to 
the "noisy nature" of AF signals. The additive noises increase the randomness of all three 
classes of signals, which makes a discriminator to judge all episodes as AF, hence favors AF 
class. As a result, the specificity for detection of AF also decreases as the SD^K reduces. 

10 This new Bayesian Discriminator has satisfactory performance (over 95%) for detection of SR, 
AFL and AF when the SNR ^ lOdB. 

To test the anti-noise performance of the disclosed discriminator, Gaussian white 
noises were intentionally added with different signal-to-noise ratio (SNR) to each episode of 
the close test data set. 

15 The effects caused by increasing the SNR on the performances of the new Bayesian 

Discriminator are presented in FIGURE 5. With a decrease m SNR, the sensitivity for 
detection of more regular rhythms as SR and AFL decreased accordingly, while the sensitivity 
for AF detection remained at high levels. However, the specificity for AF detection 
decreased with the reduction of SNR, while the specificity for SR and AFL detection remained 

20 at high levels. As a result, the overall accuracy for detection of SR, AFL and AF are similar 
at different SNRs. When the SNR is greater than 10 dB, the disclosed discriminator has an 
accuracy of about 95% in the detection of SR AFL and AF as shown in FIGURE 5. 

In addition, the presence of fiu: field R wave interfiarence also can result in 
misclassLfication of SR as AF. This problem may be addressed by, for instance, appropriate 

25 cross chamber blankiag and careful positioning of the atrial lead to avoid far field R wave may 
minimize this problem. 

5.1.7 Statistical Analysis 

Continuous variables are expressed as mean ± 1 standard deviation. The statistical 
comparisons were performed by Chi-square analysis and Student t test, as appropriate. To 
30 test the performance of the example embodiment of the disclosed discriminator, the sensitivity, 
specificity, and accuracy for detection of SR, AF, and AFL were calculated. (See Bland M. 
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In: An Introduction to Medical Statistics, Chapter 15, Oxford University Press, 1996:273-276). 
Those with P values <0.0S were considered statistically significanL 

5.1.7.1 Discriminator Performance 

The performances of the new Bayesian Discriminator for the close-test and open-test 
S data set are summarized in TABLE 2. A total of 3 episodes (4%) of &lse positive of AF 

detection occurred in 2 patients during SR due to the presence of far-field R wave sensing. 

All SO episodes of sinus tachycardia were correctly identified as SR The sensitivity, 

specificity and accuracy of riiythm detection for both close and open data set ware simiiar. 

The overall sensitivity for drtection of SR, AP and AFL is 97%, 97% and 94%, respectively; 
10 and the overall specificiQr for detection of SR, AF and AFL is 98%, 98% and 99%, 

' respectively. The overall accuracy of detection of SR, AF and AFL is 98%, 97% and 98%, 

respectively (TABLE 2). 



TABLE 2. Performances of the Bayserian Discriminator 





Rhythm Dedsioa 












SR 


AF 


AFL 


Total 


Sensitivity 


Spedficity 


Accuracy 


Close Data 


Actual rhythm 
















SR 


70 


2 


0 


72 


97.2 


98.6 


98.2 


AF 


1 


92 


1 


94 


97.9 


97.6 


97.7 


AFL 


1 


1 


51 


53 


96.2 


99.4 


98.6 


Open Data Set 


Actual iliythm 
















SR 


47 


1 


0 


48 


97.9 


97.9 


97.9 


AF 


1 


60 


1 


62 


96.8 


97.6 


97.2 


AFL 


1 


1 


33 


35 


94.3 


99.1 


97.9 
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S.L8 Main Findings 

The results demonstrate that the features of intra-cardiac atrial electrograms, which 
included the regularity, rate, en^gy distribution, percent tune of quiet interval and number of 
baseline reaching, are significantly different during SR, AFL, and AF. However, detection 

5 methods employing only one or few of these features have only limited sensitivity, specificity 
and accuracy for detection of SR, AFL, and AF. The disclosed Bayesian Discriminator 
based on the Bayes decision rule and five features of atrial electrograms, allows rapid on-line 
and accurate (98%) detection of SR, AFL, and AF with robust anti-noise performance. The 
disclosed discriminator requires a very short computing time. In an example embodiment, 

10 250ms are sufficient to make a decision for a rhythm episode of 1000ms. As shown in the 
example section, the use of multiple features discrimination provides a much higher sensitivity, 
specificity and accuracy (all >94%) for rhythm detection than single or double features 
methods, as described above. 

Clinically, as device therapies for atrial tachyanhythmias become more sophisticated 

15 in their ability to deliver several modes of therapy, such as antitachycardiac pacing and 
defibrillation, dependmg on the specific rhythm, rapid and accurate detection of potentially 
tachycardias that can be tominated by pacing will be critical Furthermore, accurate 
detection of SR fi:om AFL and AF can also prevent inappropriate device therapy. The new 
Bayesian Discriminator described m this study, which is based on nmltiple features detection, 

20 can be easy implemented in the implantable device and provides rapid (>250 msec) and 
accurate ^97%) detection of AF, with robust anti-noise performance. 

5.1.9 Conclusion 

This disclosure encompasses new methods, systems, and devices for detecting 
arrhythmias and heart diseases based on multi-variate Bayes decision, which combine a 
25 plurality of different features of the intra-atrial electrogram. The described diagnostic tools 
enable superior overall sensitivity, specificity, and accuracy for rfiythm detection than known 
single or double features methods as well as resistance to various ranges of noise. 

However, citation of documents herein is not intended as an admission that any of the 
documents cited herein is pertinent prior art, or an admission that the cited documents are 
30 considered material to the patentability of the claims of the present application. Instead, they 
are intended to clearly describe the claimed mvention. All statements as to the date or 
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representations as to the contents of these documents are based on the infonnation available to 
the applicant and does not constitute any admission as to the correctness of the dates or 
contents of these documents. 

Although the present invention has been described in considerable detail with 

S reference to certain preferred embodiments, other embodiments are possible. Therefore, the 
spirit and scope of the appended claims should not be limited to the description of the 
preferred embodiments contained herein. Modifications and variations of the invention 
described herein will be obvious to those skilled in the art fix>m the foregoing detailed 
description and such modifications and variations are intended to come withm the scope of the 

10 appended claims. Moreover, a number of references have been cited, the entire disclosures 
of which are incorporated herein by reference in their entirety. 
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What we claim: 

1 . A method for detecting a sinus rhythm of mterest with aid of a computer executable 
instructions processing device based on electrophysiological information obtained from a 
5 subject via sensors coupled to the subject, comprising the steps of: 

inputting a training data set as open-test data set of mtra-atrial electrograms for evaluation, and 
receiving a result in accordance with a set of estimated conditional probabilities of the sinus 
rhythm and a plurality of features. 

10 2. A method accordmg to claim 1 further comprising the steps of: 

selecting the phirality of features of intra-atrial electrograms and a type of output, 
inputtmg a new episode as close-test data set of intra-atrial electrograms for training, and 
estimating the set of conditional probabilities for the plurality of features and the type of 
output in accordance with a multiple-iudex Bayesian discriminator algorithm from the close- 

15 test data set. 

3 . A method according to claim 2 fiirther comprising the step of selecting additional features 
for estimating conditional probabilities. 

20 4. A method according to daim 1 fiirther comprising the step of modifyiog at least one 
estimated conditional probability firom the set of estimated conditional probabilities with the 
open-test data set and the result 

5 . A method according to claim 1 further comprising the step of diagnosing heart conditions in 
25 accordance with the result 

6. A method according to claim 2 further compridng the step of selecting a treatment to the 
subject corresponding to the result. 

30 7. A method according to claim 2, wherein the pluraUty of features is selected from the g^^ 
consisting of regularity, rate, energy distribution, percent time of quiet interval, and number of 
baseline reaching. 
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8. A method according to claim 2, wherdn the type of oirtput is selected from the group 
consisting of sinus rhythm, atrial flutter, and atrial fibrillation. 

5 9. A method according to claim 2 fiirther comprising the step of generating the result in less 
than five (5) seconds from receiving the open-test data set. 

10. A device for providing an electrical signal in response to detecting a predetermined sinus 
rhythm with thfi aid of a set of computer executable instructions based on electrophysiological 

10 information obtamed from sensors, for example, the device majLComprise: 
a set of power input terminals; 

a data module for providing conditional probabilities of the predetermined sinus rhythm 
relative to a plurality of features; 
an input for receiving electrophysiological information; 
15 a member for providing the electrical signal for a modulating effect on heartbeats; and 
at least one computer executable instructions processmg unit. 

1 1 . The device of claim 10 v^erein the predetennined smus rhythm is selected from the group 
consisting of sinus rhythm, atrial flutter, and atrial fibrillatioa 

20 

12. The device of claim 10 further comprising a module for identifying at least one feature 
from the group consisting of regularity, rate, energy distribution, percent time of quiet interval, 
and number of baseline reaching in the electrophysiological information. 

25 13 . The device of claim 10 further compriang a module nr modules for updating a conditional 
probability of the predeteoxiined sinus rhythm with the electrophysiological information. 

14. The device of claim 10 wherem the electrical signal provided by the member belongs to 
the group consisting of the electrical signal or input to the chest wall that synchronizes the 
30 heart and allows the normal rhythm to restart, small electrical impulses to the heart muscle to 
maintain a suitable heart rate, electrical energy to the heart muscle to cause the heart to beat in 
a normal rhythm, and high radio-frequency energy through a special catheter to small areas of 
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tissues that may be related to abnonnal heart rhythms. 

15 . A computer readable media carrying thereon conqniter executable instructions for carrying 
out the steps of a meUiod for detecting a sinus rhythm of interest based on electrophysiological 

5 information obtained from a subject via sensors coupled to the subject, the method comprising 
the steps of: 

inputting an open-test data set of intra-atrial electrograms for evaluation, and 

receiving a result in accordance with a set of estimated conditional probabilities of the sinus 

rhythm and a plurality of features. 

10 

16. The computer readable media of claim 15 further comprising computer executable 
instructions for carrying out the steps of: 

selecting the plurality of features of intra-atrial electrograms and a type of output, 
inputting a dose-test data set of intra-atrial electrograms for training, and 
15 estimating the set of conditional probabilities for the phirality of features and the type of 
output in accordance with a multiple-index Bayesian discriminator algorithm from the close- 
test data set. 

17. The computer readable media of claim 16 fiirther comprising computer executable 
20 instructions for carrying out the step of selecting additional features for estimating conditional 

probabilities. 

18. The computer readable media of claim 15 fiirther comprising computer executable 
instructions for carrying out the sXjsp of modifying at least one estimated conditional 

25 probability from the set of estimated conditioiud probabilities with the open-test data set and 
the result. 

19. The compute readable media of claim 15 ftirther comprising computer executable 
instracfions for carrying out the step of diagnosmg heart conditions m accordance with the 

30 result. 

20. The computer readable media of claim 16 fiirther comprising computer executable 
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instructions for carrying out the step of providing a treatment to the subject corresponding to 
the result. 

21. The computer readable media of claim 16 fiirther comprismg computer executable 
5 instructions wherein the plurality of features is selected from the group consisting of features 

such as regularity, rate, energy distribution, percent time of quiet interval, and number of 
baseline reaching. 

22. The computer readable media of claim 16 further comprismg computer executable 
10 instructions wherein the type of ou^ut is selected from the group consisting of sinus rhythm, 

atrial flutter, and atrial fibrillation. 

23. The conqmter readable media of claim 16 farther comprising conqjuter executable 
instructions for carrying out the step of generating the result preferably in less than five (5) 

1 S seconds from receiving the open-test data set 
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